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ABSTRACT 


A TOPAZ-II Thermionic Fuel Element (TFE) test stand is under study at the 
TOPAZ International Program. Extensive cross-sectional drawings of this system are 
presented along with detailed descriptions of its construction. A thermal analysis of the 
test stand 1s developed using a thermal resistance network method and data collected on 
the test stand. Critical resistances in the network include the 0.5 mm regulated helium 
gap, cooling water channel, and cesium interelectrode gap. The thermal analysis indicates 
that for an input power to the active zone of 3003 W, and a regulated helium gap pressure 
of 10 torr, the thermal conductivity of the regulated gap helium ranges from 0.073-0.087 
W/mK. This is much less than the thermal conductivity calculated as a function of 
temperature only of 0.276 W/mK. For a water mass flow rate of 0.04 kg/s, the convective 
heat transfer coefficient ranges from 528-1012 W/m*K based on correlations for fully 
developed laminar and turbulent flow. Our model may be used in conjunction with a 
thermionic emission model in order to more completely evaluate the operating 
characteristics of the test stand, and to better correlate results with actual operating 


conditions in the TOPAZ-II reactor core. 
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I. INTRODUCTION 


The TOPAZ-II Thermionic Fuel Element (TFE) test stand was originally built and 
tested at the Scientific Industrial Association "LUTCH" in Podolsk, Russia. In December, 
1992, the stand was reassembled at the New Mexico Engineering Research Institute 
(NMERI) in Albuquerque, New Mexico. The TFE test stand was designed to conduct 
single TFE non-nuclear testing by replacing the nuclear fuel with a tungsten heating 
element (TISA heater). The test stand provides a means to test a single TFE so that its 
performance can be better observed without the interference and complexity of the 
integrated TOPAZ-II system. This removes the direct coupling of the TOPAZ-II in-core 
TFE with the coolant loop, radiator and the other TFEs electrically connected in series. 
Acceptance and demonstration testing were conducted throughout 1993 (Luchau et al., 


1993). 


In order to better understand and predict the behavior of a TFE in the test stand, a 
detailed thermal analysis based on a network method needs to be developed to account for 
the thermal losses in the stand. Because of the large amounts of Russian technical data 
and engineering drawings associated with the TOPAZ-II TFE test stand, and due to the 
extreme difficulty in translating Russian technical information, the implementation of this 
thermal model required the creation of detailed and extensive cross-sectional drawings for 
this thesis. Mathcad Plus 5.0 (MathSoft, Inc., 1994) was the primary technical calculation 


software used to solve the thermal analysis model. 


Since the test stand provides a means for testing a single TFE without the support 
of the supporting TOPAZ-II reactor systems, there are several differences in the 
construction and operation of the test stand from the TOPAZ-II reactor core. In the test 


stand, water cooling is used to remove the excess heat instead of the liquid metal cooling 


that is used in the reactor core. The cooling water flows through the test stand water 
jacket along a torturous path to aid in its heat removal ability. The heat transfer 
coefficient for the cooling water is directly dependent on the flow cross-sectional area and 
the total heat transfer area, both of which are challenging to determine. Thermocouples 
internal to the test stand provide a axial temperature distribution along the TFE working 
section. No internal temperature measurements are made in the TOPAZ-II reactor core 
TFEs. Heaters are also located internally in the test stand to provide a means of degassing 
the test stand components prior to operation. A 0.5 mm gap in the test stand just inside of 
the cooling water jacket is filled with low pressure helium to provide a means of 
controlling the heat transfer through the test stand. This gap is used to help the test stand 
better simulate the operating conditions experienced by a TFE in the reactor core. Due to 
the low pressure band that the gap is adjusted over (from | to 10 torr), the pressure 


dependence of the rarefied helium's thermal conductivity is difficult to define. 


Experimental data for this work were collected in March 1994, with follow-on 
experimental data taken in August 1994 at NMERI. Along with aiding in the development 
of the thermal model for the test stand, these data are used to investigate several of the 
critical resistances from the thermal network. A good understanding of the test stand 
thermal characteristics is essential in evaluating the performance of a single TFE operating 
in the test stand, and correlating that performance to a TFE operating in the TOPAZ-II 
reactor core. The experimentally measured input power to the TFE working section along 
with the collector sleeve temperature information is also used as an input to the TOPAZ-II 
Thermionic Transient Analysis Model (TITAM). TITAM is a computer simulation 
program developed by the University of New Mexico to look at both steady-state and 
transient operating characteristics of the TOPAZ-II reactor system (Paramonov :and 


El-Genk, 1994). 


The goals of this work include the creation of detailed engineering drawings for 
the test stand, investigation of test stand material properties, development of a thermal 
model for the test stand, and benchmarking the model with the experimental data 
collected. Additionally, experimental data are used to empirically fix a range for several of 


the thermal properties that cannot be easily determined. 


LT J. R. Venable is wnting a companion thesis that will concentrate on the 
thermionic theory behind the data collected for this thesis. Venable's thesis will be 
available after March of 1995 and a joint paper will be presented ("Operational Testing 
and Thermal Modeling of a TOPAZ-II Single-Cell Thermionic Fuel Element Test Stand") 
at the Space Nuclear Power and Propulsion Symposium held in Albuquerque, New 


Mexico, on January 8-12, 1995. 





IW. BACKGROUND 


The continued exploration and exploitation of space will require increasing 
amounts of power as technology continues to advance. Supporting high definition 
television broadcasts and the rapidly growing world-wide communication needs, moving 
large level research and industrial opportunities to Earth orbit or into the cislunar region, 
revisiting and colonizing the Moon, and pushing the envelope of civilization to Mars and 
beyond are future goals which will demand large electrical power sources that are both 
compact and have long lifetimes. These larger power sources will need to survive harsh 
environments, including radiation belts, the vacuum of deep space, and hostile 
atmospheres of other planets and natural satellites in the solar system. Nuclear power can 


provide the compact power sources necessary to meet these future needs. 


A. UNITED STATES SPACE NUCLEAR POWER PROGRAMS 


Since the mid-1950s, the United States has been interested in the use of nuclear 
power in space applications. From radioisotope power systems, nuclear rockets, and the 
use of nuclear reactors for power production, interest in the space nuclear power area was 


motivated by a favorable political environment and public support until the early 1970s. 


The Transit 4A navigational satellite was launched in June of 1961 with a 
SNAP-3B Radioisotope Thermoelectric Generator (RTG) to provide power. Since that 
first successful launch of a nuclear power source, a total of 24 launches using RTGs have 
followed including the Pioneer, Viking, and Voyager missions. Although the RTG has 
proved to be very useful for many applications, the ability to support a very ambitious 
space program required the development of a nuclear reactor based on the pmnciple of 


atomic fission for both propulsion purposes and to provided electrical power. 


Project Rover was initiated at Los Alamos National Laboratory in 1955 and was 
the beginning of approximately 18 years of research and development into nuclear rockets. 
Although the nuclear rocket program was terminated in late 1972, it was the start of the 
search for a high temperature nuclear reactor. The research and demonstration test of the 
Nuclear Engine for Rocket Vehicle Application (NERVA) program proved that a nuclear 
reactor could provide a high power density rocket with high thrust and specific impulse. 


(Angelo and Buden, 1985) 


Active development of a nuclear reactor to provide electric power began in 1957 
with the Systems for Nuclear Auxiliary Power (SNAP) program. These nuclear reactors 
required a power generation method to convert the heat created by atomic fission into 
electrical power. The initial concept was SNAP-2, a thermal uranium zirconium-hydnde 
reactor that utilized a Rankine cycle with mercury as the working fluid to convert the 
nuclear heat to electricity. This was the precursor to two separate concepts, the SNAP-8 
and the SNAP-10A, both of which were initiated in 1960. The SNAP-8 reactor utilized a 
mercury Rankine cycle along with a larger version of the SNAP-2 reactor system to 


support anticipated higher power demands. (Angelo and Buden, 1985) 


The SNAP-10A was also based on the SNAP-2 reactor system, but instead of 
using the dynamic Rankine power conversion cycle, it used a thermoelectric power 
conversion system. The use of static power conversion provides a means of converting 
heat directly into electricity without the need of a secondary working fluid (such as 
mercury for a Rankine cycle or helium-xenon for a closed Brayton cycle) and the 
associated turbines and pumps/compressors. In Apnil of 1965, the SNAPSHOT mission 
was successfully launched carrying a SNAP-10A reactor to approximately a 1300°km 
orbit. Although the SNAP-10A reactor shutdown only 43 days after becoming critical, it 


proved that a nuclear reactor could be successfully launched and operated on orbit. This 


was the only attempt by the United States to launch a nuclear reactor into space. (Angelo 


and Buden, 1985) 


Another static power conversion method that was developed in parallel with the 
thermoelectric method was thermionic power conversion. Whereas thermoelectrics work 
on the basis of a thermally induced emf at the junction of two dissimilar semiconductors, 
thermionics works basically by “boiling off’ electrons from a hot emitter surface and 
collecting the electrons on a relatively cooler collector surface across a small gap. The 
collector and emitter are then electrically connected together through an external load 
which completes the electrical circuit. The process of thermionics is discussed in more 
detail in the next chapter. Although thermionic devices were tested in the mid- to 
late-1960s, no complete reactor systems using thermionics was actually built in the United 


States. (Angelo and Buden, 1985) 


Due to a changing political environment and shifting emphasis in the space 
program, all nuclear reactor development programs were canceled at the end of 1972. 
Although some research did continue after 1972, it wasn't until the early 1980s that 
increased interest was returned to space nuclear reactors. This renewed interest was 
mainly due to the need for a large power source to support the Strategic Defense Initiative 
(SDI). The advance space nuclear power program was reinitiated as the Space Power 
Advanced Reactor (SPAR) system later renamed SP-100. Several reactor designs and 
power conversion systems were initially considered for SP-100. The primary concept that 
SP-100 developed utilized a thermoelectric power conversion system. (Angelo and 
Buden, 1985) Research on SP-100 has continued until 1994 when once again the political 


environment has changed and funding for the program is being terminated. 


B. FORMER SOVIET UNION SPACE NUCLEAR POWER PROGRAMS 


The former Soviet Union, like the United States, was very active in the research 
and development of space nuclear power. From the early 1960s to 1989, a total of 38 
nuclear reactor systems were launched by the former Soviet Union. The majority of the 
nuclear reactors launched were for the Radar Ocean Reconnaissance SATellites 
(RORSATs) and utilized thermoelectric energy conversion. The other two nuclear 
reactors launched were COSMOS-1818 and COSMOS-1867. Both COSMOS satellites 
were test platforms for the TOPAZ-I reactor system utilizing multi-cell thermionic energy 


conversion. (Polansky et al., 1994) 


In the early 1960s, the former Soviet Union began development of single-cell and 
multi-cell TFEs. In 1967 the two concepts were sent for system development to two 
separate technical institutes. The single-cell TFE was developed by the Central Design 
Bureau of Machine Building (CDBMB) in St. Petersburg under the project name of 
ENISY. The multi-cell TFE was developed by the Scientific and Production Cooperation 
"Krasnaya Zvezda" in Moscow under the project name of TOPAZ, the Russian acronym 
for Thermionic Experiment with Conversion in Active Zone. Both concepts were 
developed cooperatively between the two institutes up until the early 1970s. This lead to 
many similarities between the two reactor designs. Cooperative efforts ended in the early 
1970s due to differences between the two institutes and the programs developed 


independently. (Voss, 1994) 


In 1987 and 1988, the multi-cell TFE reactor TOPAZ was successfully launched 
twice on COSMOS-1818 and COSMOS-1867. The first mission operated for 142 days 
and the second mission operated for 342 days before the planned depletion of the cesium 
supplies on board terminated the tests. Both test were extremely successful and confirmed 


that the on-orbit data correlated well to ground test results. (Paramonov and El-Genk, 


1994) Although there were 26 single-cell TFE ENISY reactor systems manufactured and 
19 systems tested between 1970 and 1989, none of the single-cell TFE technology were 


ever launched and tested on orbit. (Schmidt et al., 1994) 


C. TOPAZ INTERNATIONAL PROGRAM 


At the end of the 1980s, the Russians approached a private U.S. company 
(International Scientific Products, San Jose, CA) about transferring Russian thermionic 
technology through the purchase of Russian space reactor hardware. During the 
negotiations, the U.S. company inadvertently renamed ENISY to TOPAZ-II and TOPAZ 
to TOPAZ-I (Voss, 1994). In May of 1992, the U.S. Government purchased two 
complete TOPAZ-II reactor systems and associated test and support hardware through 
the purchasing agent International Scientific Products. Funding was provided by the 
Strategic Defense Initiative Organization (SDIO) to investigate the possible use of the 
TOPAZ-II reactor system in support of the Nuclear Electric Propulsion (NEP) Space Test 
Mission. This project has since been named the Topaz International Program (TIP). 


(Wyant, 1994) 


The TIP facility is located at NMERI which is owned and operated by the 
University of New Mexico in Albuquerque, NM. Funding is still being provided by the 
former SDIO (now the Ballistic Missile Defense Organization, BMDO) and administered 
by the U.S. Air Force's Phillips Laboratory. The New Mexico alliance for Thermionics 
which includes the University of New Mexico, Sandia National Laboratories and Los 
Alamos National Laboratory provides research support for the TIP program. The 
international effort includes scientist and engineers from the Russia, England, and France. 
Although the flight program for TIP has been canceled, in March of 1994 four additional 
TOPAZ-II units were purchased for continuing research into thermionics and reactor 


operation. (Wyant, 1994) 





Il. THERMIONIC ENERGY CONVERSION 


A. BASIC THEORY 


The thermionic process can be considered to be a "heat engine" which utilizes 
electrons as its working fluid to convert heat into electricity. The thermionic converter 
consists of a emitter, collector, and a small gap between the two. The collector is 
electrically connected to the emitter through an external load to complete the electrical 
circuit. Figure 3-1 shows a simplified flat-plate thermionic converter (Rasor, 1991). 


Electrons are "boiled" off the hot emitter and travel across the gap and are "condensed" on 
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Figure 3-1. Simplified Flat-plate Thermionic. 
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the relatively cooler collector. The potential difference between the emitter and collector 
creates a driving force which causes a current to flow through the external load. (Angelo 


and Buden, 1985) 


Since the thermionic device acts as a "heat engine," the maximum efficiency of the 
thermionic cycle can be characterized as the cycles Carnot efficiency. Some of the 
inefficiencies of the thermionic cycle include: (1) radiation heat transfer from the emitter 
to the collector; (2) negative space charge in the gap that interferes with the electron flow 
from the emitter to the collector; and (3) thermal losses. These inefficiencies prevent the 


thermionic cycle from approaching Carnot cycle efficiency. (Angelo and Buden, 1985) 


Radiation heat transfer from the emitter to the collector cannot be mitigated. 
Higher collector temperatures would not only reduced the amount of radiation heat 
transfer but would also allow for a much smaller radiator since the radiators surface area is 
inversely proportional to the radiators operating temperature to the fourth power. 
However, as the collector temperature increases, the number of electrons being "boiled" 
off the collectors surface, called electron back emission, also increases. This in effect 
causes the overall cycle efficiency to decrease dramatically. A design trade-off must be 
made to understand the optimum collector operating temperature. (Angelo and Buden, 


1985) 


The negative space charge in the gap can be decreased by either making the gap 
size smaller or by injecting a low pressure gas into the gap. A vacuum diode ts a 
thermionic device that relies on a very small gap between the electrodes to reduce the 
current inhibiting effects of the electric field produced by the negative space charge. 


Theoretical spacing needed to produce useful current density approaches the associated 
Debye length (~ | tum) which is for a practical purpose very difficult to manufacture and 
maintain. This very narrow gap can be widened somewhat (up to at most 0.5 mm) by 


injecting a gas plasma into the gap. This type of thermionic device is known as a plasma 


i 


diode. The primary purpose of this plasma is to facilitate to transfer of electrons through 
the gap while minimizing the potential drop and current attenuation due to the presence of 


the plasma. (Rasor, 1991) 


The most common gas that is used in this type of diode 1s cesium which has the 
lowest ionization potential of any chemical element (Angelo and Buden, 1985). A 
measure of the amount of energy necessary for a electron to escape through the potential 
barrier associated with the surface of a metal is called the work function and is given in 
units of electron-volts (eV) (Angrist, 1976). Not only does the cesium plasma in the gap 
decrease the space charge, it also may have the effect of forming a monolayer of cesium 
on both the emitter and collector causing a decrease in their associated work functions. 
This decrease in the metals work function is directly related to the pressure of the cesium 
in the gap and the electrode temperature. Cesium pressure can be controlled by adjusting 
the temperature of the liquid cesium reservoir which in effect controls the vapor pressure 


of the evaporating cesium. (Angelo and Buden, 1985) 


B. TOPAZ-Il NUCLEAR REACTOR 


The TOPAZ-II nuclear reactor consists of 37 individual TFEs incorporated in a 
cylindncal core arrangement. Figure 3-2 shows a cross-sectional view of the TOPAZ-II 
reactor. Each of 37 TFEs consist of a UO, fuel assemble, a single-cell cylindrical 
thermionic device, and a sodium-potassium (Nak) liquid metal cooling channel. Electrical 
power is supplied to the loads by 34 of the TFEs connected in series. These TFEs are 
designed to produce 5-6 kW, at 27 VDC for up to 3 years. The other three TFEs are 
connected in parallel to provide power to the electromagnetic pump. This pump is used to 
generate the pumping power necessary to circulate the NaK. The overall length of the 
reactor core 1s 37.5 cm with an overall width of 26.0 cm. (Paramonov and El-Genk, 


1994) 
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Figure 3-2. Cross-sectional View of the TOPAZ-II. 


The TFEs are inserted into axial channels in a zirconium-hydnde (ZrH) moderator. 
Twelve reflector drums surround the core, three are safety drums which are turned out 
prior to initial criticality when on orbit. The other nine reflector drums are used to control 
the nuclear reaction during operation and are called control drums. All of the reflector 
drums are made out of beryllium clad in stainless steel with a small section of neutron 
absorbing material made of boron-carbide (B,C). A beryllium reflector surrounds the 


remainder of the core as shown in Figure 3-2. (Schmidt et al., 1994) 


Excess heat is removed from the TFEs using a eutectic NaK coolant pumped 
through stainless steel piping by the electromagnetic pump. The hot coolant flows 
through a conical radiator where the waste heat is rejected to space. During normal 
operation the thermionic emitter temperature ranges between 1800 K to 2100 K and the 
NaK coolant maintains the collector temperature between 743 K to 843 K. (Schmidt et 
al., 1994) 
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c. TOPAZ-II THERMIONIC FUEL ELEMENT 


The TOPAZ-II TFE consists of a single-cell cylindrical thermionic converter. A 
single-cell converter utilizes the entire length of the reactor core as a single device. This 1s 
in contrast to the TOPAZ-I TFE which is made up of multi-cell thermionic converters. A 
simple comparison of the electrical circuit for both a single-cell and multi-cell TFE are 
shown in Figure 3-3. The advantages of using multi-cell converters are that they produce a 
higher power density and higher conversion efficiency. This is mainly due to less ohmic 
losses in the multi-cell because of the shorter thermionic elements. Another advantage of 
the multi-cell design is that each of the cells are connected in series resulting in a TFE 
Output voltage in the volt range. This reduces the amount of electrical losses in the leads 


and makes the power distribution system more efficient. (Paramonov and El-Genk, 1994) 
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Figure 3-3. Simplified Electrical Schematics for Multi-Cell and Single-Cell TFEs. 
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The two major disadvantages of using a multi-cell TFE are that manufacturing is 
significantly more complicated and the design does not provide a means of removing the 
nuclear fuel and performing non-nuclear testing. Although the single-cell TFE design has 
a much lower efficiency, it does provide a less complicated converter and a means of 
performing non-nuclear testing. This testing is accomplished by replacing the UO, fuel 
assembly with a tungsten heater (TISA heater). (Paramonov and El-Genk, 1994) The 
TISA heaters provide a temperature distribution similar to that created using nuclear 
heating. Detailed information on the use of electrical heaters to simulate nuclear heating 


can be found in El-Genk and Xue, 1994. 


The interelectrode gap between the emitter and collector is filled with cesium 
vapor during normal operation. The cesium is supplied by a single liquid cesium reservoir 
that is heated by one of the NaK coolant return lines. The cesium in the gap ts constantly 
vented to space at a rate of approximately 0.5 g/day. This is one of the limiting factors in 
the reactor lifetime. (Schmidt et al., 1994) During normal operation, the cesium pressure 


in the gap is maintained around 1.5-2.0 torr (Paramonov and El-Genk, 1994). 


Axial and radial temperature distributions in the reactor core cause performance 
differences in the 37 TFEs. The exact performance characteristic of each individual TFE 
is difficult to determine due to the following: (1) the thermal coupling of each TFE with 
the TFEs surrounding it, or lack of TFEs; (2) the electrical coupling of each TFE with the 
other TFEs, either in series for the 34 load power producing TFEs or in parallel for the 3 
pump power producing TFEs; (3) all TFEs share a common heat removal system; and (4) 
all TFEs share a common cesium system. To understand the individual thermionic 
performance of a TFE, the TOPAZ-II TFE test stand must be used. This test stand allows 


detailed analysis of a single TFE under unique operating characteristics. 


IV. TEST STAND DESCRIPTION 


The TOPAZ-II TFE test stand is used to perform non-nuclear testing of a single 
TOPAZ-II TFE. This allows the TFE performance to be determined without the 
interference of the rest of the TOPAZ-II reactor systems and without any interfering 
effects caused by the TFEs relative location in the core (see Figure 3-2). Interference by 
the reactor systems include the direct thermal coupling of the TFE to the NaK coolant, 
coolant pump, and radiator, and the electrical coupling of the TFE to the other TFEs 
connected in series in the reactor core. The test stand also allows varying certain 
parameters which could not be varied on the entire reactor system to investigate their 
effect on TFE performance. These parameters include varying cesium vapor pressure in 
the interelectrode gap, the removal of all cesium vapor from the interelectrode gap, 
changing the collector sleeve temperature (directly effecting the collectors temperature), 


and operating the TFE at its operational limits to determine break-down characteristics. 


The test stand contains the test apparatus within a vacuum chamber along with the 
associated vacuum, cesium, helium, and cooling water systems external to the vacuum 
chamber. As with almost all Russian structural components, the vacuum chamber is 
constructed mainly of Russian stainless steel type 12X18H10T. Phillips Laboratory has 
determined that this type of Russian stainless steel is closely approximated by the 
American Iron and Steel Institute (AISI) standard 321. The vacuum system is used to 
initially evacuate the vacuum chamber, the helium cavities, the cesium vapor system, and 
the TISA heater channel. The test stand start-up procedure requires an involved 
evacuation and out-gassing procedure. Out-gassing is performed by numerous heaters 


throughout the test stand. 
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The cesium system 1s used to supply cesium vapor to the interelectrode gap and to 
maintain the cesium vapor pressure. The system also provides a means to purify the 
cesium inventory. When the test stand is operating, cesium vapor is provided to the 
interelectrode gap by a cesium reservoir that is inside the vacuum chamber. This isolates 
the internal cesium system from the rest of the system outside of the vacuum chamber. 
Cesium vapor pressure in the interelectrode gap is regulated and varied by adjusting the 
inner cesium reservoir saturation temperature. An external cesium reservoir 1s provided to 
evacuate and fill the cesium system when the test stand is not in operation. (Luchau et al., 


1993) 


The helium system provides helium to the cavities and gaps inside the test stand, a 
means to purify the helium supply and a means to evacuated the helium from the system. 
Helium is used to fill in the gaps in the test apparatus to facilitate and control the heat 
transfer process. Collector sleeve temperatures can be controlled by varying the pressure 
in the regulate helium gap. The mechanism by which the thermal conductivity of helium is 


affected by pressure is discussed in detailed in the next chapter. 


The water cooling system picks up the excess heat throughout the test stand. The 
components cooled by the water system include electric drives for various valves, vacuum 
pumps, cesium reservoirs, transistor load bank and the test apparatus. The volumetric 
flowrate of the cooling water flow to the test apparatus is measured using a turbine type 
flowmeter in the inlet line. The inlet and outlet cooling water temperature is measured by 
type-K, surface contact thermocouples on the inlet and outlet stainless steel piping. Two 
thermocouples are used on the outlet piping, one on the outlet pipe that is in thermal 
contact with the vacuum chamber and the other on the outlet pipe that is not in thetmal 
contact with the vacuum chamber. The experimental data show that both of these 


temperature measurements are well within the thermocouple tolerance of each other. The 
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flowrate and temperature measurement accuracies are critical to the thermal modeling of 


the test stand. 


The stand utilizes an automatic control system, an automated data acquisition 
system, and a variable electrical load bank. The entire control and data acquisition system 
was completely upgraded from the original Russian system to a United States technology 
based system. This upgrade occurred between the first set of tests in March of 1994 and 
the second set in August of 1994; however, the testing parameters were maintained the 
same during both data sets. The load bank 1s composed of both a transistive and a 
resistive load bank. The load bank can be controlled by a micro-computer to provide 
current-voltage (I-V) sweeps for the TFE. These sweeps are important in determining the 


actual operating charactenistics of the TFE. (Luchau et al., 1993) 


A diagram of the test apparatus is shown in Figure 4-1 (based on Luchau et al., 
1993). The upper and lower helium chambers provide helium to the unregulated gaps in 
the test section and to the spaces in the TISA heater. The cesium for the inter-electrode 
gap enters the test section through the lower helium chamber which is heated during 
testing to prevent "cold-spots" in the cesium line. Since the lower helium chamber 1s 
heated and the upper helium chamber is not heated, heat losses out of the top and bottom 
of the test apparatus are different. This difference in the end losses will affect the test 
stand's thermal model. The upper and lower helium chambers also contain the electrical 
connections to both the emitter and collector. Electrical connections to the TISA heater 


are made in a chamber above the upper helium chamber. 


A lengthwise cross-sectional view of the test stand's working section is shown in 
Figure 4-2. This figure shows the relative proportions of the various test section 
concentric cylinders. The overall length of the thermionic working section is 375.0 mm. 


The length of the TISA heater is 300.0 mm and the length of the cooling water jacket 1s 
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Figure 4-1. Cross-sectional View of the TOPAZ-IJ TFE Test Stand. 
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approximately 326.0 mm. The extremely small dimensions and close tolerances gives an 
indication of the difficulty associated with the manufacturing of the test section. A radial 
cut-away view of the complete test stand is shown in Figure 4-3. An enlarged view of the 
portion of the test stand from the collector sleeve inward is shown in Figure 4-4. The 
overall diameter of the entire test section is only 87.0 mm giving an L/D of 4.31 for the 


working section. 


The TFE and its associated TISA heater are fit into the collector sleeve tube of the 
test stand as shown in Figure 4-4. The TISA heater is constructed of a special tungsten 
heating element and flexible lobes that act as thermocompensators. The maximum supply 
voltage 1s 29 VAC with a maximum current of 170 A. The TISA heater can supply a 
maximum heating power to the working section of 4500 W (Wold et al., 1993). The inner 
conductor of the TISA heater has a diameter of 6.5 mm and the outer conductor has a 
thickness of 0.4 mm and an outer diameter of 7.0 mm. The length of the heater is 300.0 
mm which makes the heater 37.5 mm shorter at each end than the actual thermionic 
working section. This shorter length simulates the end effects associated with nuclear 
heating. Detailed discussions of the correspondance of TISA heating to nuclear fuel 


heating are presented in El-Genk and Xue (1993). 


The thermionic device consists of concentric emitter and collector cylinders 
separated by an interelectrode gap and is 375.0 mm in length. The emitter cylinder is 1.15 
mm thick with an outer diameter of 19.6 mm. The emitter is made of a monocrystalline 
molybdenum - 3% niobium alloy with its exterior surface coated with a 0.1 mm thick 
monocrystalline tungsten layer. The 0.5 mm interelectrode gap is filled with cesium vapor 
during normal operation. The interelectrode gap is maintained by the use of 
scandium-oxide (Sc,O,) spacers which prevent shorting between the emitter and collector. 
A total of six spacer sets are located axially along the TFE working section. These 


spacers are not shown in any of the figures. The collector cylinder is 1.4 mm thick with an 
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Figure 4-4. Radial Cut-away Detail of a TOPAZ-II TFE in the Test Stand. 
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inner diameter of 20.6 mm. The collector is made of polycrystalline molybdenum with its 
outer surface coated bya 0.15 mm layer of alumina (AI,O,). The alumina electrically 
insulates the TFE from the remainder of the test mg. The thermionic device and TISA 
heater used in the test stand are identical in construction to that used in the actual 


TOPAZ.-II reactor. 


A 0.05 mm gap separates the alumina insulator from the collector sleeve tube. 
This gap allows the TFE to be fit into the test ng and is filled with helium during normal 
operation. The helium is provided by the upper and lower helium chambers (see Figure 
4-1). The helium pressure in this gap is unregulated and is normally in the 200 to 300 torr 


range. 


The collector sleeve tube is 3.0 mm thick with an outer diameter of 29.9 mm. It is 
constructed of Russian 1X18HI1OT stainless steel. The collector sleeve contains 12 evenly 
spaced, 2.0 mm deep slots in its outer surface that run the length of the collector sleeve. 
Figure 4-5 shows the relative location of the thermocouple slots in the collector sleeve. 
These slots contain the only thermocouples located in the interior of the test ng. The 
thermocouples are positioned along the collector sleeve to provide an axial temperature 
distribution in the working section. The thermocouple wires are packed in alumina inside 
of stainless steel tubes which are laid inside the slots. These thermocouples are positioned 
in the test ng at the location where the NaK liquid metal cooling flows through the TFE in 
the actual TOPAZ-II reactor. This allows correlation of the test stand data to data 
obtained from the complete reactor testing. Another 0.05 mm helium gap separates the 
collector sleeve tube from the copper conductor. This small gap is for assembly purposes 
only. The helium in this gap is also supplied from the upper and lower helium chambers 


and is unregulated in pressure. 
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The copper conductor is 15.0 mm thick with an outside diameter of 60.0 mm. It is 
made of 99.99% pure copper and provides a heat transfer medium outward to the cooling 
water. The outer surface contains five radially spaced slots that each contain four 
nichrome (NiCr) heating elements. These heaters are used for outgassing the rig prior to 
startup and are secured during normal operation. The heater wires are insulated with 
alumina and sheathed in stainless steel. A 0.5 mm gap separates the copper conductor 
from the inner cooling water jacket. This gap ts filled with helium supplied from a source 
external to the test mg. The pressure of the helium ts regulated in the band of | to 10 torr 
during normal operation to control the collector sleeve temperature by varying the thermal 
resistance of the gap. The exact mechanism through which the helium thermal 


conductivity is affected by helium pressure is discussed later in the thesis. 


The cooling water jacket 1s made of Russian 1X18HIOT stainless steel and has an 
inner wall thickness of 2.5 mm and an outer wall thickness of 1 mm. A coiled stainless 
steel tube between the jacket walls provides a spiral flow path for the cooling water. The 
water flows from the bottom to the top of the water jacket around the outside of the coil 
in the cooling "channel" depicted in Figure 4-2. The coils are spaced by 35.0 mm and 
spiral approximately 6.5 times from bottom to top. Cooling water temperature is 
measured on the inlet and outlet piping exterior to the vacuum chamber using surface 
thermocouples. Cooling water flow rate is measured in the inlet piping using a 


turbine-type flowmeter located exterior to the vacuum chamber. 
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V. TEST STAND THERMAL ANALYSIS 


As discussed in Chapter III, the thermionic cycle is a relatively inefficient 
conversion process. The TOPAZ-II TFE has a maximum efficiency of 5-6%; therefore, at 
least 94% of the heat generated by the TISA heater must be dissipated as a loss. Detailed 
thermal analysis of the test stand can provide a better understanding of how testing a 
single TFE in the test stand correlates to operating a TFE in the TOPAZ-II reactor core. 
The developed thermal model also yields valuable information about the limitations of the 


test stand for future testing and the validity of past and future test results. 


A. THERMAL MODEL INPUTS 


The amount of input power to the test stand is determined by measuring the input 
current to the TISA heater through an external shunt and by measuring the voltage drop 
across the input terminals of the TISA heater. Multiplying the current in amps by the 
voltage in volts gives the input power to the TISA heater in watts and is designated as 
Pisa. his TISA input power must be corrected to calculate the amount of power actually 
delivered to the emitter of the TFE. Correction is required to take into account the lead 
and line losses associated with the heater. Although Russian documentation for 
calculating the correction factor was not available, the Russian technicians at TIP use a 
correction factor of 88%. Correcting the input power to the TISA heater by this factor 
provides the amount of actual power delivered to the "active zone" of the TFE. This 


active zone power is designated as P,.. 


The output power produced by the TFE is determined by measuring the current 
through the external load bank and by measuring the voltage drop across the upper and 


lower emitter or collector terminal connections (see Figure 4-1). Since the voltage drop is 


Zo 


measured directly across the interior terminal connections of either the emitter or 
collector, no correction factor is required to determine the output voltage of the TFE. 
Multiplying the voltage in volts by the current in amps gives the TFE output power in 


watts and is designated as P.... 


The only temperatures that are actually measured internally in the test stand are the 
collector sleeve temperatures. These measurements occur 1.8 mm into the stainless steel 
collector sleeve from its inner cylindrical surface (see Figure 4-4). As discussed in 
Chapter IV, the thermocouples in the collector sleeve provide axial temperature 
measurements along the length of the working section. There are a total of 12 
thermocouples in the collector sleeve which are designated T56 through T67. During the 
data collection runs, thermocouple T64 was inoperative and was disregarded in the 
analysis. Thermocouples T59 and T60 and thermocouples T61 and T62 measure the 
temperature at the same location along the collector sleeve. Relative locations of the 
thermocouples with respect to the centerline of the working section are shown in Figure 
5-1. The lengths of the TISA heater, the thermionic working section, and the water jacket 
and the approximate locations of the cooling water inlet and outlet are also shown in 


Figure 5-1. 


In a one-dimensional model, it is necessary to take an average of the 11 operating 
collector sleeve thermocouples to obtain a value to use as an input temperature. Although 
many variables effect the magnitude of the collector sleeve axial temperature profile, the 
profile shape is generally the same under all conditions. Figure 5-2 shows a typical axial 
temperature profile at a TISA heater input power of 2400 W. The temperature data 
points are shown as diamonds with T56 at the -205 mm position and T67 at the 205'mm 


position. The dashed line represents a third-order polynomial curve fit of the data. 
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Figure 5-2. Typical Axial Temperature Distribution along the Collector Sleeve. 


Although the TISA heater provides a fairly uniform temperature profile into the emitter, 
several reasons inherent to the test stand design cause the distribution to take the shape as 
shown in Figure 5-2. An overall downward slope in the data is expected due to the 
cooling water flowing from the bottom to the top of the working section; however the end 
effects of the upper and lower helium chambers cause the data to not be linear at the two 
ends of the graph. A straight average of the 11 data points would bias the temperature 
input to the thermal model low due mainly to the end effects. This biasing can be 
minimized if the end data points (T56 and T67) are not used and an average is taken only 
of the remaining nine data points. For example, the average of the 11 data points in 
Figure 5-2 is 379.0 K. Excluding the end data points gives an average of 392.1 K which is 
approximately 3.5% higher. For all calculations performed in this analysis, the average 
collector sleeve temperature was determined by excluding T56 and T67. A much more 
detailed two-dimensional model would be required to understand the error associated with 


this approximation. 


a2 


The cooling water system provides the last of the measured test stand quantities 
necessary for the thermal analysis. Instrumentation of the cooling water system was 
required in support of this thesis since the information on the cooling water system was 
not previously deemed necessary. Surface thermocouples were placed on the cooling 
water inlet and outlet lines, exterior to the vacuum chamber. Approximately one meter of 
piping separates the actual thermocouple location from the inlet and outlet piping 
penetrations on the test stand. No corrected factor was used on these temperature 
measurements for any potential piping losses between the test section and the exterior of 
the vacuum chamber because the only method of heat transfer from the water piping 
would be by radiation heat transfer to the interior of the vacuum chamber. Since the 
maximum water temperature inside the piping is less than 40 °C over all the data runs, 
radiation heat transfer is not significant. The water piping is also in thermal contact with 
the vacuum chamber which provides another potential heat loss by conduction heat 
transfer to the vacuum chamber itself. Some uncertainty could have been removed if the 
thermocouples had been placed on the piping inside the vacuum chamber at the entrance 
and exit of the test section. This was not done for any of the data runs used in this 
analysis due to time constraints and the complexity of removing the vacuum chamber. An 
average bulk water temperature in the water jacket can be estimated by taking an average 


of the inlet and outlet water temperatures. 


A turbine-type flow meter was also installed on the inlet cooling water line prior to 
the data runs. The accuracy of the flow meter was physically verified to be within 2% of 
the indicated flow rate over the range used during the data runs. Initially the data runs 
were taken with the cooling water isolation values fully open which was the normal 
operating procedures for the test stand. However, with the valves fully open, the change 


in cooling water temperature between the inlet and the outlet was less than 5 °C. Since 
the accuracy of the surface thermocouples is only + 1 °C, a larger temperature change is 
desired to reduce the percent error. Around halfway through the first set of data runs, the 


flow rate was throttled down to approximately 0.7 gpm. This gave a temperature change 
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between 10 to 15 °C over most of the data runs. All of the data collected during the 


second set of data runs were collected at a flow rate of approximately 0.7 gpm. 


B. THERMAL MODEL 


Difficulty was encountered during the thermal analysis of the test stand due to the 
small number of input parameters and the large number of "unknowns." In order to better 
understand the scope of the problem, a thermal analysis flowchart was developed. Figure 
5-3 shows the test stand thermal analysis flowchart. Data collected from the test stand are 
shown as circle and are used as inputs into the thermal model. From this flowchart it can 
be seen that the problem divides itself into the following three sections: (1) thermal 
network calculations; (2) water channel heat removal calculations; and (3) coupling the 
thermal network calculations with the TOPAZ-II Thermionic Transient Analysis Model 
(TITAM). TITAM its a computer model developed by the University of New Mexico to 
simulate the operation of the entire TOPAZ-II reactor system under steady-state and 
transient conditions (Paramonov and El-Genk, 1994). The basis for the entire analysis is 


the thermal network. 


1. Thermal Network 


The one-dimensional losses through the test stand are modeled using a network 
method. Figure 5-4 shows a simplified representation of the resistance network used for 
thermal modeling the test stand. The network starts at the outside wall of the cylindrical 
emitter and runs outward to the centerline of the water channel. The three major modes 
of heat transfer encountered throughout the network are radiation, conduction, and forced 
convection. Radiation heat transfer was considered a heat transfer mode through all of the 
test stand gaps when the thermal network was initially developed. It is shown later in this 


section that radiation heat transfer plays an insignificant role in the heat transferred 
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through the unregulated and regulated helium gaps, however, radiation heat transfer does 
play an important role in the cesium interelectrode gap. Conduction heat transfer occurs 
through all of the metal components and through the gaps in the test stand; including the 
cesium interelectrode gap, the two unregulated helium gaps, and the regulated helium gap. 


Forced convection occurs into the cooling water flowing through the water channel. 


The basic heat transfer equation used for conduction heat transfer through a 
cylindrical wall is found in Wong (1977). Equation 5.1 gives the equation used for the 
metal components of the test stand: 


Q a (T hot = lizoia) (s 1) 


1 
R. 


where: R. = thermal resistance (K/W) 
Tyo = hot wall temperature (K) 
T .oig = Cold wall temperature (K) 


The thermal resistance, R., is calculated for hollow cylinders using Equation 5.2: 


In(=2 


= Ss 
2nzk 2) 





where: r, = outer wall radius (m) 
r, = inner wall radius (m) 
z = cylinder length (m) 
k = material thermal conductivity (W/mK) 


For the test stand, T,,, will be the temperature at the inner cylindrical surface and T,,,, will 
be the temperature at the outer cylindrical surface of the component of interest. For the 
thermal model calculations, the temperature subscript will indicate the node at which the 


temperature 1s measured. 


oy) 


A cylindrical length of 30.0 cm 1s used in all of the calculations. This length 
corresponds to the total length of the TISA heater; however, as can be seen in Figure 5-1, 
the length of the TISA heater, the length of the water jacket, and the length of the TFE 
working section are all different. Since there is a 20% increase in cylinder length from the 


TISA heater length to the TFE working section length, the thermal resistance also 
increases by 20%. Although this only results in a few degree Kelvin difference in the AT 


across a single metal component, the cumulative effect of this difference can lead to large 
calculational uncertainty. The shortest of the possible cylinder lengths is used in the 


analysis since it gives a more conservative result. 


The thermal conductivities of all the different test stand components are graphed as 
a function of temperature in Appendix A. Since the thermal conductivity is a function of 
temperature and the metal components have a decreasing temperature distribution from 
the inner wall to the outer wall, an estimate of the thermal conductivity for the entire 
thickness must be made. When entering the thermal model, only one temperature is 
known; therefore, it is not possible to determine the average temperature of the 
component. For the first iteration through the network, the thermal conductivity of the 


component is calculated based on the one known temperature. Because metal is a very 
good conductor of heat, the thermal conductivity is fairly constant over the small AT 
effective across the metal. Based on calculations, using the one known temperature to 


estimate the thermal conductivity results in only a 0.4% difference than if the thermal 
conductivity were recalculated with the lower temperature found after calculating the AT. 
Therefore, less than a three degree Kelvin error is induced into the network and this makes 


recalculation of the thermal conductivity unnecessary. 
The basic equation used for heat transfer for the gas filled gaps is found in El-Genk 


et al. (1993). This equation accounts for both radiation and conduction heat transfer 


through the gap and Is given in Equation 5.3: 
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Q = 7 (T hot wi T cold (5.3) 


r/¢ 


where: R.. = radiation/conduction thermal resistance (K/W) 


The radiation/conduction thermal resistance, R,., 1s calculated using Equation 5.4: 


To 
In(y) Ke 2 2 =| 
{ dad SE f(T hort T cot (Thot+T cota) (5.4) 


al 202(T> — Ti) 


where: k, = gas thermal conductivity (W/mK) 
d = gap width (m) = (r, - 1.) 
o = Stefan-Boltzmann constant = 5.67x10° W/m°K"* 
Eer = effective emissivity 


The thermal conductivity for helium as a function of temperature is graphed in Appendix 
A. According to Reid et al. (1977), the thermal conductivity of all gases can also be 
affected by the pressure in the gap, although for pressures greater than approximately 1 
torr, the thermal conductivity only increases by about 1% for a pressure change of 760 
torr. Therefore, for the unregulated helium gap pressures of 200-300 torr, the thermal 
conductivity is considered a function of temperature only. This estimation is valid because 
the mean free path for collisions between atoms 1s much less than the mean free path for 
collision between an atom and the walls; therefore, an atom does not have to travel far 
before transferring its energy to another atom and eventually reaching thermal energies. 
The pressure maintained in the regulated helium gap is around the lower end of where this 
approximation is considered valid. According to the experimental data collected, the 
thermal conductivity of helium in the 1-10 torr range is greatly affected by pressure. This 


unusual phenomenon Is discussed in detail later in this chapter. 


The effective emissivity can be calculated by assuming that the view factors from 
both surfaces are equal to 1, and the two surfaces areas are also equal by Equation 5.5 


(Welty et al., 1984): 
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oe ae (5.5) 


where: €, = inner surface emissivity 
€, = outer surface emissivity 


Due to the extremely small dimensions and gap size, both of these assumptions are valid. 


The basic equation used for forced convection heat transfer into the cooling water 
is found in Wong (1977). The equation used for calculating the heat transfer into the test 


stand cooling water is based on Newton's Law of Cooling and is given in Equation 5.6: 
eae 
Q= Ra (Twait ~ Toutk) (5.6) 


where: R,. = forced convection thermal resistance (K/W) 
T ya) = Cooling water channel wall temperature (K) 
Ty. = bulk cooling water temperature (K) 


The forced convection thermal resistance, R,, is calculated using Equation 5.7: 


In 


a Piyed (c= 75 oe 


Re 
where: h = convective heat transfer coefficient (W/m’K) 


The convective heat transfer coefficient is a complicated function of the cooling channel 
geometry, the thermophysical properties of the cooling water, and principally of the flow 
characteristics through the channel. This coefficient and the complexities encountered in 


its definition are discussed in detail later in this chapter. 
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An example of the data set used in the thermal network calculations is given in 
Table 5-1. These data are used in the following discussion to explain the thermal model 
calculational procedure. The data collected between 14-18 August, 1994, in support of 
this thesis are tabulated in Appendix B. The data presented in Table 5-1 were taken on 
August 17, 1994. As can be seen there, the TFE efficiency is less than 4%. This is mainly 
due to the interelectrode cesium pressure not being at its optimum value for an input 
power of 3000 W. At this power level, an interelectrode cesium pressure of 
approximately | torr would increase the TFE's efficiency closer to 5%. However, for 
thermal modeling purposes, optimum TFE performance is not important. Cesium 
optimization will be covered in detail in the companion thesis wntten by LT J. R. Venable 


that will be available after March of 1995. 





Ps, (W) 3,412 Minimum Collector 354.3 
Sleeve Temp. (C) 







P., (W) 3,003 Maximum Collector 5156 
Sleeve Temp. (C) 
P4@y) 109 Average Collector 453.7 


Sleeve Temp. (C) T, 


Q,., (W) 2,894 Cooling Water Flow 0.69 

rate (gpm) 
Efficiency (%) 3.63 Cooling Water Mass 0.04 

Flow rate (kg/s) 
Cs Pressure 0.5 Cooling Water Inlet 16.6 
(torr) Temp. (C) | 
Regulated He 10 Cooling Water Outlet 32.2/32.1 |. 
Gap Press (torr) Temp. (C) 


Table 5-1. Example Data Set for the Thermal Model. 
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Since the only temperature measured internally to the test stand 1s the collector 
sleeve temperature, the thermal network is entered at node T,. From this node, 
calculations can be made to determine the temperature distribution either inward toward 
the center of the TFE or outward toward the water channel. In order to make these 
temperature calculations using basic heat transfer equations, the amount of heat 
transferred out through the test stand must be calculated. This transferred heat is equal to 
the portion of the TISA input power that is not converted to electricity in the TFE and ts 


designated as Q,..,. Equation 5.8 shows the definition for Q....: 


Qhoss == Pw tg P out (5.8) 


In a one-dimensional model for the test stand, Q,.. will travel only radially outward 
through the test stand components to be carried away by the cooling water flowing 


through the water jacket. 


Using the node 6 temperature (T,) and the amount of heat lost by the system (Q....) 
as a Starting point, the thermal model is presented in two sections: (1) calculations from 
the collector sleeve to the water channel, and (2) calculations from the collector sleeve to 
the emitter. Detailed discussions about the critical resistances in the network are also be 


presented. The calculations made with the data from Table 5-1 are given in Appendix C. 
a. From the Collector Sleeve to the Water Channel 
The first resistance encountered in the thermal network is a 1.2 mm thick 
stainless steel collector sleeve. The thermal conductivity of the stainless steel ts 


determined using Figure A-3. This thermal conductivity is then used to calculate the 


thermal resistance of the stainless steel using Equation 5.2. Using the calculated Q,,., as 
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the variable Q, Equation 5.1 can now be solved for T_,,, which corresponds to T,, since all 


of the other equation variables have been defined. 


A 0.05 mm unregulated helium gap is the next resistance in the thermal 
network. The thermal conductivity of this continuum helium gas (gap pressure of 
approximately 200 torr) is calculated using Figure A-4. The effective emissivity is then 
calculated using Equation 5.5. This equation requires the stainless steel emissivity which 
is found from Figure A-6 and the copper emissivity which is found from Figure A-7. Both 
of these emissivities are fairly constant over the temperature range of interest. Using 
Equation 5.4 in combination with Equation 5.3, the value for T, can be solved for 
numerically. If the contribution of the radiation heat transfer is ignored through the helium 
gap, and T, is solved for based on conduction heat transfer through the gap only, the 
answers for T, are the same within four significant digits. For this reason, radiation heat 
transfer through the unregulated and regulated helium gaps can be disregarded and the 


resistance network can be solved based on conduction heat transfer only. 


The next resistance in the network is the copper conductor. Again 
Equations 5.2 and 5.1 are solved to determine the temperature on the outer cylindrical 
surface of the copper conductor, which corresponds to T,. The regulated helium gap 1s 


the next resistance in the thermal network. 
b. Regulated Helium Gap 


As discussed earlier in this chapter, the pressure of the helium in the 
regulated gap is maintained between 1-10 torr, depending on the input power level, to 
hold the collector sleeve temperature fairly constant. Since the collector sleeve 
temperatures are measured at approximately the same radial location that the liquid metal 


cooling flows through the actual TOPAZ-II core, holding this temperature constant better 
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simulates the actual conditions in the core. The temperature at this radial position in the 
TOPAZ-II core is relatively constant due to the high heat capacity of the liquid metal 


cooling and the high heat rejection capability of the radiator. 


Due to the low gap pressure, the helium is rarefied which indicates that the 
mean free path for collisions between atoms is much larger than the mean free path for 
atom to wall collisions. Paul Agnew of AEA Technology, United Kingdom, currently 
assigned to the TOPAZ-II Project and working in the area of materials science and 
thermionics, has completed a preliminary analysis on the pressure dependence of helium 
thermal conductivity at these low pressures. The initial results of his work indicated that 
helium thermal conductivity should approach its continuum value at a gap pressure of 
approximately 5 torr vice the 1 torr discussed in Reid et al. (1977). However, 
experimental results indicate that the helium thermal conductivity should be around 30% 
of its continuum value for a gap pressure of 10 torr. Detailed Monte Carlo computer 
programming and investigation of gas atom scattering mechanisms will be required to 


determine a more accurate model of the regulated gap. (Agnew, 1994) 


Although a more exact solution is not possible without extensive work, the 
regulated gap helium thermal conductivity can be better approximated using the empirical 
data collected on the test stand. Since the regulated helium gap is separated from the 
cooling water channel by only 2.0 mm of stainless steel, the temperature of the inner wall 
of the cooling channel is directly dependent on the regulated helium gap thermal 
conductivity. Using the helium continuum thermal conductivity value calculated to be 
0.276 W/mK for this example, gives a cooling channel inner wall temperature of around 
325 °C. The cooling water is pressurized on average to a pressure not greater than 70 
psig. Since the cooling system provides water to many different components, the actual 
pressure 1s dependent on the cooling system load conditions. This cooling water pressure 


gives a initiation boiling point temperature of approximately 140 °C. The wall temperature 
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calculated using the helium continuum thermal conductivity is around 185 °C greater than 
this boiling point temperature. Based on an cooling water outlet temperature of 32 °C, it 
is assumed that the water does not boil in the cooling channel. Figure 5-5 shows a graph 


of the cooling water channel inner wall temperature dependence on the regulated helium 
gap thermal conductivity. If the wall temperature is bracketed around 80 + 30 °C, the 
regulated helium gap thermal conductivity ranges between 0.073-0.087 W/mK. This 


range is around 30% of the calculated continuum value. 


This calculated range of helium gap thermal conductivities can be used 
along with Equation 5.1 to calculated a range of node 10 temperatures. Using these node 
10 temperatures, the 2.0 mm thick stainless steel water jacket thermal conductivity can be 


determined using Figure A-3. The temperature at the inner wall of the cooling water 
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Figure 5-5. Cooling Water Channel Inner Wall Temperature Dependence on 
the Regulated Helium Gap Thermal Conductivity. 
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channel, corresponding to the node 11 temperature, is then calculated using Equations 5.2 


and 5.1. These calculated wall temperatures range from 64-122 °C. 


Using an estimated convective heat transfer coefficient and an estimated 
heat transfer area (the method used to estimate these parameters is discussed later in this 
chapter), the bulk water temperature can be estimated using Equations 5.7 and 5.6. These 
calculations yield a range of bulk water temperatures of 6-92 °C. Using an average of the 
measured inlet and outlet cooling water temperatures gives a bulk temperature of 


approximately 24 °C. 
c. From the Collector Sleeve to the Emitter 


The calculations for the inward portion of the thermal network are similar 
to the calculations made for the outward portion. The first resistance encountered is the 
inner portion of the stainless steel collector sleeve. The same thermal conductivity value is 
used for the collector sleeve as was used in the outward portion. Using Equations 5.2 and 


5.1, the node 5 temperature can be found. This corresponds to T,. 


The next resistance in the network is another 0.05 mm unregulated helium , 
gap. Radiation heat transfer is ignored and the heat transfer calculations using Equations 
5.2 and 5.1 provide the temperature at node 4. Node 4 corresponds to the outer wall of 
the layer of alumina insulation on the collector. The thermal conductivity of alumina 1s 
found using Figure A-land the same conduction heat transfer calculations are completed 
to calculate the node 3 temperature, T,. This temperature is the interface temperature 


between the collector and the alumina insulation. 


The temperature at node 2, T,, can be similarly calculated using the 


thermal conductivity of polycrystalline molybdenum found from Figure A-5. This 
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temperature corresponds to the inner surface of the collector and is an extremely 
important measurement to determine thermionic performance. The cesium interelectrode 


gap 1s the next resistance in the network. 
d. Cesium Interelectrode Gap 


The cesium interelctrode gap heat transfer mechanism is another critical 
resistance in the network that cannot be easily calculated with only basic heat transfer 
equations. Conduction and radiation heat transfer both play an important part in the 
overall heat transfer through the interelectrode gap. The thermal conductivity of the 
cesium in the gap depends on the gases vapor pressure, the emitter surface temperature, 


the collector surface temperature, and the gap size. (Paramonov and El-Genk, 1994) 


Radiation heat transfer cannot be ignored through the interelectrode gap 
since the temperature of the emitter 1s high enough to make this mode of heat transfer 
significant (on the order of 1800 K). The effective emissivity for the gap cannot be 
calculated as presented for the helium gaps, but must be determined empirically. A 
detailed discussion on thermal modeling of the cesium interelectrode gap and coupling the 
thermal model with a thermionic emission model can be found in Paramonov and El-Genk 


(1994). 
2. Water Channel Calculations 


In order to calculate the thermal resistance using Equation 5.7, the convective heat 
transfer coefficient is required. Calculating the heat transfer coefficient for the water 
channel 1s complicated by the flow channel geometry and by the low water flow rates. 
Figure 5-6 shows the cooling water jacket details including the flow guidance coil. This 


flow guidance coil is a stainless steel tube with an outer diameter of 10.0 mm that fits 
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Figure 5-6. Cooling Water Jacket Details. 
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inside of the water jacket to direct the flow in a spiral path to facilitate forced convection 
heat transfer. The cooling water flows around the outside of the flow guidance coil with 


no flow through the coil itself. 


As shown in Figure 5-6, the flow channel cross-sectional area is calculated by 
determining the cross-sectional area of a rectangle with dimensions of 35x10 mm, and 
subtracting off the cross-sectional area of the flow guidance coil (diameter of 10 mm). 
This gives a cross-sectional channel flow area of 2.7x107° m*. The equivalent hydraulic 


diameter is calculated using the definition from Wong (1977), to be 0.011 m. 


Following the definitions and equations presented by El-Genk and Rao (1986), the 
Reynolds number ts calculated using this geometry and the given example flow rate to be 
1480. This Reynolds number indicates that the flow falls into the combined laminar flow 
regime. The convective heat transfer coefficient can be defined from El-Genk and Rao 


(1986) and is given in Equation 5.9: 





De (5.9) 


where: Nu = Nusselt number 
k = water thermal conductivity (W/mK) 
D. = equivalent diameter (m) 


In the combined laminar regime, the Nusselt number is dependent on the thermophysical 
properties (defined by the Graetz number) and the buoyant forces (defined by the Rayleigh 


number). This flow regime is immensely difficult to analyze with existing correlations. 


Although the convective heat transfer coefficient cannot be determined exactly, 
several existing forced convection correlations can provide a range of possible 


coefficients. The upper limit is determined using the Dittus-Boelter correlation for fully 
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developed turbulent flow in a channel. This gives a convective heat transfer coefficient of 
1012 W/mK. The lower limit is found by the Seider-Tate correlation for developing 
laminar flow. This correlation gives a convective heat transfer coefficient of 528 W/m°K. 
Therefore, the actual coefficient will fall between the range of 528-1012 W/m’K and these 
limits can then be used in the resistance network to estimate the drop between the cooling 


channel wall temperature and the bulk water temperature. 


Using the cooling water inlet and outlet temperatures along with the measured 
cooling water flow rate, the amount of heat carried away by the water can be calculated 


using Equation 5.10: 


Q =m Cp (Tout a Ta) (5.10) 


where: m = cooling water mass flow rate (kg/s) 
c, = specific heat capacity of water 
To = Cooling water outlet temperature (K) 
T,, = cooling water inlet temperature (K) 


The heat carned away by the water, Q, should be equal to the heat lost by the TFE, Q,..., if 
all of the excess heat traveled one-dimensionally out of the test section to be carried away 
by the cooling water. For the example data, the heat carried away by the cooling water is 
approximately 2840 W. This is less than 2% different from the 2894 W lost by the TFE. 
As can be seen from the experimental data presented in Appendix B, most of the data 
collected for the heat removed by the cooling water falls within 10% of the heat lost by 
the TFE. The percent error tends to increase as the input power level decreases. This 1s 
an expected result since the difference between the inlet and outlet cooling water 
temperature decreases as the power level decreases causing an increase in the expected 


CIror. 
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By setting Equation 5.10 equal to Equation 5.6, the convective heat transfer 
coefficient can be calculated based on a heat transfer surtace area. Once again difficulty 
arises in determining the exact heat transfer surface area. From Figure 5-6 it can be seen 
that even though no cooling water flows through the flow guidance coil, the stainless steel 
coil is in thermal contact with the stainless steel water jacket. The maximum possible heat 
transfer surface area can be determined by summing the inner water jacket cylindrical 
surface area, the outer water jacket cylindrical surface area, and the total surface area of 
the flow guidance coil. The minimum possible heat transfer surface area is equal to the 
cylindrical surface area of the inner water jacket. These calculations are shown in 
Appendix D and give a range of heat transfer area between 0.067-0.203 m*. Figure 5-7 


shows a graph of the heat transfer coefficient plotted against the heat transfer area. Since 
T(wall) = 45 °C 
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Figure 5-7. Cooling Water Heat Transfer Coefficient Dependence on the 
Heat Transfer Area. 
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the cooling water bulk temperature is known, the cooling channel wall temperature 1s 
varied from 45 °C to 95 °C. It can be seen from the figure that a wall temperature of 
greater than 95 °C would fall outside of the heat transfer coefficient/heat transfer area 
bands. This information can be used along with Figure 5-5 to further limit the regulated 


helium gap thermal conductivity used in the thermal resistance network. 


3. Thermal Network Coupling With TITAM 


The TITAM computer code provides a means of simulating the operation of a 
TOPAZ-II TFE in the reactor core under both steady-state and transient conditions 
(Paramonov and El-Genk, 1994). Although the test stand is constructed differently than 
the actual core and the thermal and electrical coupling of the TFE with the actual reactor 
systems does not exist, the location of the collector sleeve thermocouples in the test stand 
does provide a means of entering TITAM and having TITAM determine the operating 


characteristics of the TFE. 


The required inputs into the electrically heated TFE component operation portion 
of TITAM are thermal input power, inlet NaK coolant temperature, and NaK coolant flow 
rate. The thermal input power to the TFE is equal to the active zone power, P,,, used in 
the thermal network. The average collector sleeve temperature, T,, may be used as the 
NaK coolant inlet temperature, and a typical NaK flow rate of 1.3 kg/s is used as the last 
required input to TITAM. From the inputted data, TITAM outputs the radial temperature 
distribution through a TFE in the reactor core and provides the output operating 
characteristics. The TITAM program was executed using the inputs from Table 5-1. The 
steady-state output and the relevant data from the thermal analysis are presented in Table 


5-2. 
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Parameter TITAM Output | Thermal Analysis 
Results 
Thermal Power to 3,003 3,003 
TFE (W) 
(V) 
Load Current 184.38 78.89 
(A) 
Load Electric Power 149.58 108.87 
(W) 
TFE Efficeincy 4.98 3.63 
(%) 
Cesium Reservior 576 528 
Temperature (K) 
Emitter Surface 1835557 
Temperature (K) 
Collector Surface 765.71 763.93 
Temperature (K) 
Insulator WL a be 761.88 
Temperature (K) 
Coolant Average 727.58 726.83 
Temperature (K) 


Table 5-2. Comparison of TITAM output to Thermal Analysis Results. 











The discrepancy in the TFE output power could be the result of the differences 
between operating a TFE in the test stand and simulating the operation of a TFE in the 
reactor core. Another possible source of error are the initial and operating conditions that 
TITAM uses in its calculations being different from the actual test stand conditions. 
TITAM has a user fmendly interface through two input files which can be changed to 
reflect various initial and operating conditions in the reactor core. This includes the ability 


to vary the steady-state interelectrode gap cesium pressure. 
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TITAM initially uses a steady-state cesium pressure of approximately 2.0 torr 
(based on a cesium reservoir temperature of 576 K). The input files were not altered prior 
to running the example data set through TITAM to obtain the above results. When 
TITAMs initial steady-state cesium pressure was changed to reflect the actual operating 
pressure in the test stand of 0.5 torr (based on a cesium reservoir temperature of 528 K), 
the TFE output power increased from 151 W to 188 W. Although TITAM is being used 
for a purpose that it was not designed for, this trend is contrary to the experimental data 


collected. 


The TITAM temperature distribution correlates well with the thermal analysis 
mainly due to these calculations being based on simple heat transfer principles. This 
indicates that TITAM could be used to estimate the emitter surface temperature based on 


test stand thermal analysis data. 
In order to understand the discrepancies encountered by trying to couple the 


experimental test stand data with TITAM, a detailed examination of how TITAM actually 


performs it calculations and the assumptions under which it operates is required. 
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VI. CONCLUSIONS 


The TOPAZ-II TFE test stand provides an effective means to investigate the 
performance of an individual TFE without the interfering effects and complexity of the 
integrated TOPAZ-II system. A detailed thermal analysis is completed to better predict 
the behavior of a single-cell TFE and to understand the thermal losses associated with the 


test stand components. 


Prior to developing the thermal model, extensive research of the test stand 
construction and materials was necessary. Due to the large volume of Russian 
documentation and the difficulty in translating highly technical data and engineering 
drawings, the development of test stand cross-sectional drawings and the determination of 


component materials and material properties were crucial to the thermal analysis. 


A thermal resistance network has been utilized to determine the thermal loss 
mechanisms through each of the test stand components. Conductive heat transfer through 
the metal components accounts for a very small portion of the overall temperature drop 
from the emitter to the water channel. For the data set used in the example calculations, 
only 7% of the estimated temperature drop from the inner surface of the collector sleeve 


to the bulk water temperature occurred in the metal components. 


Both radiation and conduction heat transfer occur through the regulated and 
unregulated helium gaps, but calculations indicate that radiation heat transfer does not 
contribute significantly to the total heat transferred through these gaps. The pressure in 


the regulated helium gap is varied during power changes to maintain the collector sleeve 
temperature relatively constant. Due to the very low pressure maintained in this gap (< 10 
torr), the helium is rarefied and the thermal conductivity of the gas is not only dependent 
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on gas temperature but is also a function of gas pressure. Estimations based on the 
assumption that the cooling water is not boiling in the channel give a range of regulated 
helium thermal conductivities between 0.073-0.087 W/mK at a gap pressure of 10 torr. 
This range is approximately 30% of the thermal conductivity value calculated assuming a 
temperature dependence only. The regulated helium gap is one of the most critical 


resistances in the thermal network. 


The primary mode of heat transfer into the water channel is forced convection heat 
transfer. Determining the amount of heat transferred into the water by forced convection 
requires calculating the convective heat transfer coefficient. This coefficient 1s dependent 
on the flow characteristics in the channel, the cooling channel geometry and the 
thermophysical properties of the cooling water. The Reynolds number is calculated to be 
1480 which indicates that the flow rate is in the combined laminar flow regime. This 1s a 
very difficult regime in which to accurately calculate the heat transfer coefficient due to 
the importance of the buoyant, conductive and convective forces in the channel. Based on 
existing correlations, the convective heat transfer coefficient is determined to fall between 


the range of 528-1012 W/m’K. 


Another critical resistance in the thermal network is the cesium interelectrode gap. 
Conduction heat transfer occurs through the cesium vapor and can be defined in terms of 
the cesium vapor pressure, the gap size and the surface temperatures of the emitter and 
collector. Due to the high emitter temperatures, radiation heat transfer is significant 
through this gap. The gap heat transfer characteristics are directly affected by the 
thermionic process. Coupling the thermal model with a thermionic model is crucial to 


estimating emitter temperatures based on data collected from the test stand. 


The TITAM computer code provides a means of simulating steady-state and 


transient operation of a TOPAZ-II TFE in the reactor core. TITAM can be utilized to 
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estimate the operating characteristics of a TFE in the test stand by inputting the active 
zone power, average collector sleeve temperature (corresponding to the Nak 
temperature), and an assumed Nak flow rate (based on typical values). The collector 
sleeve temperature calculated by TITAM is 766 K compared to a temperature of 764 K 
calculated using the thermal network. These temperatures calculations are based on 
simple heat transfer principles which provides good correlation between the two 
calculations. TITAM can also be used to estimate the emitter surface temperature based 
on the input data; however, the thermal resistance of the cesium interelectrode gap must 
be defined before the thermal analysis presented in this thesis can be used to provide an 


estimate of the emitter surface temperature to compare it to. 


Differences in the TITAM output and the thermal analysis results can be attnbuted 
to several reasons. The most important of these is that TITAM simulates a TFE operating 
in the TOPAZ-II reactor core and not in the test stand. Several other differences include 
cesium reservoir temperature, TFE load voltage, TFE load current, and TFE efficiency. 
To understand the discrepancies between the two results requires a detailed examination 
of the TITAM code and integration of the thermal model developed in this thesis with the 


thermionic model of the code. 
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VII. RECOMMENDATIONS 


The TOPAZ-II TFE test stand provides an excellent vehicle for investigating the 
thermionic characteristics of a single TFE; however, understanding the thermal 
characteristics of the test stand and how its thermal behavior affects the thermionic 
process is crucial in being able to correlate test stand results to the entire reactor core. 
Several improvements to the test stand and additional testing could greatly aid in 


improving the thermal model developed in this thesis. 


A detail investigation of the flow characteristics of the cooling water would 
provide a better estimation of the convective heat transfer coefficient. This could include 
moving the inlet and outlet cooling water thermocouples inside the vacuum chamber 
closer to the test sections water jacket. Using more accurate temperature measurements 
would allow increasing the water flow rate into a better understood and documented flow 
regime. Testing of the cooling water characteristics could be accomplished using a known 
power input into the degassing heaters and would not require the operation of the TFE. 


Thermocouple instrumentation inside the water channel itself, would also prove useful. 


Further understanding the pressure dependence of the regulated helium gap's 
thermal conductivity is critical to the accuracy of the thermal model. Once the thermal 
characteristics of the regulated helium gap and the cooling water channel are well defined, 
a more accurate heat balance can be performed between the collector sleeve and the 


cooling water. 


Exploring the end effects at the upper and lower helium chambers would provide a 
better understanding of the 2-dimensionality of the thermal system. Experiments could 
easily be performed without the operation of the TFE to examine how these end chambers 


are affecting the axial temperature distribution along the collector sleeve. Operating the 
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heaters for one or both of these chambers during operation of the TFE provides an 
additional heat input into the thermal model; therefore, the input power to these heaters 


should be recorded. 


Experimental information is only valuable if the accuracies of the instrumentation 
and data acquisition system are completely defined. The integration of Russian and 
American equipment is inherent to the TOPAZ project and increases the difficulty in 
developing an accurate model for the test stand. Improved documentation of the test 


stand measurement errors should accompany any future experimental work. 


Experimental data should be collected from the test stand at a predetermined TFE 
output voltage. This voltage level should be based on supplying a specified bus voltage 
(i.e. 28 VDC or 30 VDC). This must be a consideration since operating the TFE at the 
maximum output power (or maximum efficiency) might not produce useful voltage to 


supply a power distnbution bus. 


TITAM is a valuable tool that should be modified further for use with the test 
stand. Utilizing the thermionic emission model from TITAM along with the thermal 
model developed in this thesis can provide a detailed understanding of the thermal 
characteristics of the test stand. This would allow for a more direct correlation of the test 


stand results with the operation of the overall TOPAZ-II reactor system. 
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APPENDIX A. THERMAL PROPERTIES 


The thermal conductivities and emissivities of various test stand components are 
graphed and are defined in terms of temperature. Most of the information was available in 
Paramonov and El-Genk, 1994, but some was derived from reference book data tables. 
References are provided in the titles of each figure. A polynomial curve fit of the 


tabulated data accomplished using Mathcad Plus 5.0 (MathSoft, Inc., 1994). 
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Figure A-1. Thermal Conductivity of Alumina (AlI,O,). 
(Paranomov and El-Genk, 1994) 
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Figure A-3. Thermal Conductivity of AISI-321 Stainless Steel. 
(Peckner and Bernstein, 1977) 
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Figure A-5. Thermal Conductivity of Polycrystalline Molybdenum. 
(Paranomov and El-Genk, 1994) 
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Figure A-7. Emussity of 99.99% pure Copper. 
(Baumeister, 1978) 


1000 





1100 


68 


OOTl 





(8L6] ‘Jojslouneg) 
‘(fo'Ty) eurumypy Jo Assia “g-y amsty 


(3{) sINyeJodurs |, 
0001 006 008 OOL 009 


(L- S-Erst)-, Wy O1-r961 + 870= (Lys 


00S 


OOP 


b0 


SPV 0 


¢ 0 


gS 0 


D0 


69 





APPENDIX B. EXPERIMENTAL DATA TABLES 


Data were collected in support of this thesis during two separate experimental 
sessions. The first occurred between 15-26 March, 1994, when data were collected over a 
TISA heater input power range of 1000 W to 3600 W. The primary purpose of the runs 
completed during this period was to establish the operating characteristics of the test stand 
and its associated TFE from normal operating range to very low input power levels. This 
information is being used in a follow-on thesis to examine the thermionic characteristic of 


a TFE in the test stand. 


The second experimental session occurred between 14-18 August, 1994. Pnor to 
this session, the entire test stand data acquisition system was converted from a computer 
hardware and software system that was furnished by the Russians to a more contemporary 
system that was designed and assembled by the US, TIP personnel. This system update 
provides for much easier data acquisition and test stand control. Although the data 
acquisition system changed between the two experimental procedures, the test stand 
Operating conditions were maintained as close as possible to allow for data correlation 
between the two collection periods. The pnmary purpose of the second experimental 
procedure was to support thermal model calculations by investigating the relationship 
between regulated helium gap pressure and collector sleeve temperature and to collect 
data a various power levels for an optimized cooling water flow rate. Relevant data to the 


thermal analysis are tabulated in this appendix. 


The appendix is separated into three time periods, each time period consisting of 
three tables. The shaded columns contain data that were recorded from the test stand and 
the unshaded columns contain derived data. All of the data are collected at a cesium 


reservoir temperature of approximately 255 °C (around 0.5 torr of pressure). 
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The second column of the first table is TISA active zone power and 1s determined 
by multiplying the TISA input power by 88%. The fifth column is TFE output power and 
is determined by multiplying the TFE output voltage by its output current. The column 
indicated as column "A" is the power loss from the TFE. These values are found by 
subtracting the TFE output power from the TISA active zone power. TFE efficiency in 
column seven is the result of dividing the TFE output power by the TISA active zone 


power and multiplying by 100. 


Column eight of the first table is the average collector sleeve temperature and 1s 
the same value as in column five of the third table. Performing a heat balance on the 
cooling water provides the data in the columns indicated as column "B" and "C". There is 
one column for each of the two output thermocouples. These values are the same values 
as given in columns 12 and 13 of the third table. The % error in columns 11 and 13 of the 
first table are the percent difference between the cooling water heat balance (columns "B" 


and "C") and the power balance (column "A") 


The second table presents the 11 collector sleeve thermocouple measurements. 
See Figure 5-1 for the relative axial location of the thermocouples along the collector 
sleeve. Columns 2 through 6 of the third table provide relevant information concerning 
the thermocouple measurements. Column five gives the average of all 11 collector sleeve 
thermocouples and column six gives the average of thermocouples T57 through T66. This 


excludes the end point temperature measurements of T56 and T67. 
Column eight of the third table contains the measured water flow rate (in gpm) 


converted to a cooling water mass flow rate (in kg/s). Columns 12 and 13 of the third 


table indicate the results of the cooling water heat balance (utilizing Equation 5.10). 


WZ 


The tables dated 13-14 August 1994 are taken over a TISA input power of 
1000-3600 W. At each power level the regulated hellum gap pressure was adjusted to its 
optimum value. This maintains a relatively constant collector sleeve temperature over the 
entire power range. The high % error for the 1000 W data run is due to the small 
temperature difference between the cooling water inlet and outlet. This indicates how 


sensitive the cooling water heat balance 1s to this temperature difference. 


The tables dated 14-16 August 1994 provides data collected over three power 
ranges with varying regulated helium gap pressure at each power level. The first table 
shows that for a constant input power level, the collector sleeve temperature decreases as 
the regulated helium gap pressure increases. This proves how important understanding 
the pressure dependence of the regulated helium gap thermal conductivity is to the thermal 


analysis of the test stand. 
The tables dated 17-18 August 1994 provides a set of data runs under similar 


conditions as presented in the first set of tables. A complete set of all data acquired during 


both experimental periods is archived at the TIP facility in Albuquerque, New Mexico. 
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APPENDIX C. THERMAL NETWORK CALCULATIONS 


The following calculations and the calculations in Appendix D were completed using 
Mathcad Plus 5.0 (Mathsoft, Inc. 1994). This calculational software uses a colon followed 
by an equal sign (:=) to mean "defined by" vice an equal sign (=) to mean "equal to". All 

variables must be defined prior to being used in an equation. After a variable is defined, an 
equal sign can be used to determine the variables value. 


THERMAL MODEL INPUTS 
1. INPUT POWER TO THE TISA HEATER 
Q TISA ° 3412-watt 
2. POWER TO THE ACTIVE ZONE 
Qar = QqTs4:0.88  Q,, =3.003-10° -watt 
3. TFE OUTPUT POWER 
4. POWER LOSS FROM TFE 
3 
Qioss = Qaz~ Pout RQloss = 2-894°10 watt 
5. RADIAL DIMENSIONS FOR THE TEST STAND 
ry: 0.0098-m rg = 0.0119-m rz = 0.01495-m 
r> = 0.0103-m r5 = 0.01195-m_ reg = 0.0150-m 


rz = O001175-m  re¢ = 0.01370-m = rg = 0.0295-m 


6. LENGTH OF THE TEST SECTION 


z = 0.3-m 


7. DIMENSIONAL CONSTANT FOR CALCULATIONS 


Co 2-1-7 c = 1.885-m 
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r 79 = 0.0300-m 
ry] = 0.0325-m 


ry = 0.0375-m 


8. COLLECTOR SLEEVE AXIAL TEMPERATURE DISTRIBUTION 
T 56 = 424.49K Tso - 483.9-K Te = 433.73-K Teg = 362.95-K 


T 57 = 515.64-K = Teg = 482.57K  T63 = 439.32-K Tea = 354.29-K 
T 5g = 49846K Tg, - 490.12K Ts = 376.42-K 
9. AVERAGE COLLECTOR SLEEVE TEMPERATURE 


The average collector sleeve temperature is calculated without using the end point 
temperatures, T(56) and T(67). 


157+ T5g+T59+ Teo + 161 + 162+ 163+ 165+ T 66 


Teer 
6 9 


C 


Tg = T6c+273.15K  T6 =726.829-K 


B. FROM THE COLLECTOR SLEEVE OUTWARD 


1. STAINLESS STEEL COLLECTOR SLEEVE 


a. Thermal Conductivity of Stainless Steel 


watt 


watt att 
k = 12.0-;-——} + 0.01. 
SS bed 


WwW 
Jie kK gg = 19.268 -— — 
m-K? K 











b. Temperature Drop Across Collector Sleeve 


[7 
Q xs 
'6 





c. Temperature at node 7 
ls, ‘= age = AT gg Ty = Ne Sez he 
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2. UNREGULATED HELIUM GAP 


a. Thermal Conductivity of a Helium Continuum 











/ \ 
k Her > 9376107)". >- 913.107 Sp,’ 
mK] \m:K°, 
tt -13 | watt | 
ee AC a Ty - 1.26410. |——)-T 3" 
\m-K m:K’ , 








-2 | watt \ - watt 
tee OG ewig ieee teem Ose 
m:-K m-K 


b. Temperature Drop Across the Unregulated Helium Gap 


ig 


Q ‘In 
loss : 


/ 
AT He 7 a AT Hea 18.147-K 
He 





c. Temperature at Node 8 


T¢ > 1 > Pale Tg = /01.725-K 
3. COPPER CONDUCTOR 


a. Thermal Conductivity of Copper 


watt 


watt 


m-K? 





igaiy = 422.286, oe 6 71409.10 
m:K 
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b. Temperature Drop Across the Copper Conductor 


it g\ 
Q An| — | 
loss \rg| 


ck Cy 








AT cy = AT ¢, = 2.767-K 


c. Temperature at Node 9 


T9 - Tg aa en T 9 = 698.958-K 


4. REGULATED HELIUM GAP 


The thermal conductivity of the helium in the regulated gap is first estimated using 
the equation for a continuum gas. 


a. Thermal Conductivity of a Helium Continuum 



































-4 | watt -7 | watt 2 
kHe3 + 9376-10 °-|———|-Tg - 9.1310 | ———|-T 9 
m:-K m:-K 
-10 / watt 3 -13 | watt 4 
ae [ene ie) ae g 
m:K m:K 
-2 | watt _ watt 
K Hereg = -9.06-10 ("3 qe kK He3 as kK He4 K Hereg = 0.210 sae 


b. Temperature Drop Across the Regulated Helium Gap 


r10 


9 


Q Joss In 








Wgee = AT Hereg = 93:366-K 


CK Hereg 
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c. Temperature at Node 10 Based on a Helium Continuum in the Gap 
Tiocont = Pg - AT Hereg T 10cont = ©09.592-K 


This temperature will lead to a water channel inner wall temperature of around 
325 °C. Since this is approximately 200 °C higher than the boiling point of water, the 
helium thermal conductivity must be much lower than its continuum value. From Figure 
5-5, the helium thermal conductivity is estimated to be 0.08 W/mK for a wall temperature 
of 80 °C. Using a temperature range of 80 + 30 °C gives a low thermal conductivity of 
0.073 W/mK and a high thermal conductivity of 0.087 W/mK. The temperature drop 
across the regulated helium gap 1s recalculated based on the low, estimated, and high 
thermal conductivities. 


d. Estimated Temperature Drop Across the Regulated Helium Gap 





10 
Fas 
k O07s ee = : a 
Helow «0 KAT Helow * TAT Helow 7 353-441-K 
°K Helow 
r'10 
Q loss In ne 
watt 
KHeest = 0-08—— AT peep = ———— + AT pest = 322-515-K 
m-K C-K Heest 
ake 
Q toss nl —— 
k 0.087 at aT AT = 296.565-K 
enh = <-VO% Hehigh” Hehigh “7 > = 
mK ¢’K Hehigh 


e. Estimated Temperature at Node 10 


Tiolow = Fg - AT Helow qT 10low = 345.517-K 
T 10 = T9 7 AT Heest fh 10 = 376.443-K 
T 1Ohigh * Tg AT Hehigh T 10high = 402.393-K 
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S. Stainless Steel Water Jacket 


a. Thermal Conductivity of Stainless Steel 














watt watt _ watt 
K wilow = 120) + 001 5} 10low k wjlow = 15.455 aed 
m:K™ | 
tt 
ee 12.0 001g =a k yj = 15.764. 
m:K mike m:K 
watt watt watt 





b. Temperature Drop Across the Water Jacket 











ri 
Q joss n r10 
wjlow 
al 
Q Joss’ !| — 
AT = ‘10 AT = 7.795-K 
veal ok. oy 
W) 
ri] 
Q Joss In 
AT .. ct: ee AT |. -1:-1 = 7.668-K 
high high 
wyhig Ckenen wyhg 


c. Temperature at Node 11 


Tillow = T 1olow ~ AT wjlow T | low = 337.567-K 


Tythigh = Tiohigh- AT wjhigh — T 1 thigh = 394.724-K 
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These calculations yield a water channel inner wall temperature of between 
64-122 °C. Using the Dittus-Boelter correlation convective heat transfer coefficient of 
1012 W/m°’K as the upper limit and the Seider-Tate correlation convective heat transfer 
coefficient of 528 W/m°’K as the lower limit, the last resistance in the network can be 
calculated and a bulk cooling water temperature can be calculated. 


6. WATER CHANNEL 


a. Convective Heat Transfer Coefficient and Heat Transfer Area 








2 
trae - 528- = Aatlow - 0.07-m 
| 2 
mK A thigh = 0.125-m 
watt 
hhich = 1012: Pehla: fonhich 
: m=-K A ht : Aaa a thig A ht = 0.098-m? 


Z 


b. Temperature Drop from the Wall to the Bulk Fluid Temperature 








Q loss 

AT a ee ae AT slow 56.209-K 
ht low 
Q loss 

AT whigh ° rea" : AT whigh = 29.327-K 
ht" high 


c. Bulk Water Temperature 


T bulklow1 ad Lllow ~ AT wlow T pulklowl = 281.357-K 


T bulklow2 = ! lllow ~ 4 whigh T pulklow2 = 308.24-K 


Denier =) i Al view T pulk1 = 312.439-K 
T pulk? - T co AT whigh T Rul = 339.322-K 


T pulkhighl = Tithigh- 4? wiow — T bulkhight = 338-515-K 
Toulkhigh2 ~ Tithigh- AT whigh 7 bulkhigh2 = 365.398-K 


These calculations give a bulk cooling water temperature range of 8-92 °C. 
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FROM THE COLLECTOR SLEEVE INWARD 


1. STAINLESS STEEL COLLECTOR SLEEVE 


a. Thermal Conductivity of Stainless Steel 


watt 


watt 
Te — kgg= 19.268. 


m-K? 


t 
ban 12.0“ + 0.01. 











b. Temperature Drop Across Collector Sleeve 


6 
Q| wi 
OSS \t 5 


ok ss 








AT cg: = AT oo; = 10.888-K 


c. Temperature at node 5 


Ts = T6+ AT co: T 5 =737.717-K 


2. UNREGULATED HELIUM GAP 


a. Thermal Conductivity of a Helium Continuum 


























-4 | Watt -7 | watt 2 
sae - 9376-10 - ; Ts = oO lsaom- : oi 
m-K m-K 
; tt ’ 
k Heaj = 5.5-10°-|——|-T 5° - 1.264-10°"* = Ts 
m-K m:-K 














-2 / watt 7 watt 
k pe; = - 9.06-10 ("3) + kpei + kK yer; k Hej = 0-288 _— 
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b. Temperature Drop Across the Unregulated Helium Gap 


won 
Ctoss a | 


4/ 
© K Hei 
c. Temperature at Node 8 


T4-T5+ATye Tq =760.099-K 


3. ALUMINA INSULATION 


a. Thermal Conductivity of Stainless Steel 


Pa Gee RT, SS Oe 
| 4 4 
re 2 3 


m:K m:K 


e 2 





alum 


-§ watt 3 
en? - -7,.78-10 - pla + 1.31-10 


m:K m-K? 


-]] Watt 4 
—— Ty 





210916 oe 
m:K 


k alum ~ k aluml ~ k alum2 k alum 


b. Temperature Drop Across the Alumina Insulation 


'4 
Q os 
AT oo AT 


= —— = 1.784-K 
c-k 


alum ~ alum 


alum 
c. Temperature at node 3 


3 - Tg + AT T 3 = 761.883-K 


alum 
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4. COLLECTOR 


a. Thermal Conductivity of Polycrystalline Molybdenum 





watt -2 watt watt 
k - 110-—— - 1.5-10 ~.——.T k = 98.572 -— — 
polyMo mine Bae 5 polyMo ag 
b. Temperature Drop Across the Collector 
3 
Q Joss I 2 
a ees 7 AT eo = 2-051-K 
©°K nolyMo 
c. Temperature at Node 2 
T> > 13 Al oy T 4 = 763.934-K 
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APPENDIX D. WATER CHANNEL CALCULATIONS 


WATER CHANNEL HEAT TRANSFER AREA 


1. Dimensional Constants 


Tepe eco ee Inner Cylinder Radius in m 

fe OO425 my =o Outer Cylinder Radius in m 

z=-0326-m —  ——— qnee- Water Jacket Length in m 

ee SOOO S]My eee Radius of the Flow Guidance Coil 

eee OFO Tm fares Test Stand Diameter to center of the Water 
Channel in m 

8 - 13-deg = =  ——— ---2- Cant Angle of Flow Guidance Coil in degrees 

N=-65 eee Number of Revolutions of the Flow Guidance 
Coil 


2. Inner Cylinder Surface Area in m?: 
a = 2-it-t 11°2 

A; = 0.067-m" 

3. Outer Cylinder Surface Area in m?: 


Ag = 21 13-2 





ed 2 
A Aa 0.087-m 
4. Surface Area of the Flow Guidance Coil in m?: 

dq 
d an = Si Effective Test Stand Diameter based on the Flow 
cos(@) Guidance Coil Cant Angle in m 

d app = 0.077-m 
L-Nwdop nse Length of the Flow Guidance Coil in m 
L =1.572-m | | 
ee 4 0ek , ----- Circumference of the Flow Guidance Coil in m 
C , =0.031-m 
Ags LC. we Surface Area of the Flow Guidance Coil in m? 
A. =0.049-m° 
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5. Total Surface Area in Water Channel in m?: 


eur eie(a | fe TaN yee JAX 


A, =0.203-m" 


HEAT TRANSFER CORRELATIONS 


1. Definitions of Physical Properties of Water at 15°C (Welty et al., 1984): 





k 
o = 1000.0- = ----- Density 
m 
joule . 
c, = 4185.0-;——} a= Specific Heat Capacit 
: bees ; _— 
k a ° ° ° 
u = 11,6-10 ‘(| Dynamic Viscosity 
m-sec 
2 a 
yen nGaT one m_ ----- Kinematic Viscosity 
sec 
t bags 
k = oss, *) ----- Thermal Conductivity 
m-K 
Cyl 
Pr = . Pr =8.228 _ ----- Prandtl Number 


2. Total Cooling Water Flow Area in m’: 
A ¢ = (0.035-m)-(0.010-m) - =-(0.005-m)° 
Ae=2.715-10° +m 

3. Wetted Perimeter in m: 
Pw = (2:0.035-m) + (2-m-r 6] 
P, =0.101-m 


oe 


4. Equivalent Hydraulic Diameter in m: 


ai 
Dee ee 
Pw 
Dae oO Ollem 
5. Cooling Water Flow Velocity in m/s: 
ae 0,69.| 84 | V =4353-10° «msec. === Measured Volumetric 
\mun/ Flow rate 
V 
v=— 
aN 


v =0.16-m-sec | 


6. Reynolds Number: 





D,v'p 
Re - : 
Ul 

Re = 1.48-10° 


7. Seider-Tate Relation for Laminar Flow: 


Rew = Keb Pe=1.218-10) ----- Peclet Number 
x6 kg : ; : Br: 
eS 52:10 ee eee Wall Temperature Fluid Viscosity for 80°C: 
m-sec 
3 0.14 
k D : 
h ST = * 16(re 5] > 
be L Lw/ 
watt 





h oT = 528.245 - 5 
m'-K 


8. Dittus-Boelter Correlation for Fully Developed Turbulent Flow: 





k 
Minne == -0.023-Re Pr 
De 
3 watt 
h pp =1.012-10° - 
pp =1.012-10° -— 
m=-K 
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